Experimental evidence for the breakdown of a Hartree-Fock approach in a weakly 

interacting Bose gas 
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We investigate the physics underlying the presence of a quasi-condensate in a nearly one dimen- 
sional, weakly interacting trapped atomic Bose gas. We show that a Hartree Fock (mean-field) 
approach fails to predict the existence of the quasi-condensate in the center of the cloud: the quasi- 
condensate is generated by interaction-induced correlations between atoms and not by a saturation 
of the excited states. Numerical calculations based on Bogoliubov theory give an estimate of the 
cross-over density in agreement with experimental results. 
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Since the first observation of the Bose-Einstein conden- 
sation in dilute atomic gases, self consistent mean-field 
approaches (see e.g. 0,12) have successfully described 
most experimental results. A good example is the criti- 
cal temperature at which the condensate appears [1,0,13 ■ 
However, the approximate treatment made in these the- 
ories about particle correlations can not capture all the 
subtle aspects of the many body problem and the success 
of these theories relies on the weakness of interactions 
in dilute atomic gases {pa^ ^ 1 where p is the atomic 
density and a the scattering length for three-dimensional 
(3D) gases) and also on the absence of large fluctuations 
(Ginzburg criterion) Q. In a 3D Bose gas, this latter cri- 
terion is not fulfilled for temperatures very close to the 
critical temperature Q and deviations from the mean- 
field theories are expected. For a given atomic density, 
the critical temperature is calculated to be slightly shifted 
towards a larg er value as compared to the mean-field pre- 
diction 0, lid. ITU but this discrepancy is still beyond the 
precision of current cold atom experiments 0, 0| . 

In mean-field theory, a Bose gas above condensation is 
described by a Hartree-Fock (HF) single fluid in which 
atomic interactions are taken into account only via a 
mean-field potential V^.f. = 2,gp0,Q where g is the cou- 
pling constant and p the atomic density. Interaction- 
induced atomic correlations are neglected and the gas is 
modeled as a group of non-interacting bosons that expe- 
rience the self-consistent potential Vm.i.- As for an ideal 
Bose gas, the two particle correlation function at zero 
distance g^^HO) is 2 (bunching effect) This HF single 
fluid description holds until the excited state population 
saturates, which is the onset of Bose condensation. 

In this paper, we present measurements of density pro- 
files of a degenerate Bose gas in a situation where the 
trap is very elongated and the temperature of the cloud 
is close to the transverse ground state energy. For suffi- 
ciently low temperatures and high densities, we observe 
the presence of a quasi-condensate 151 at the center of 
the cloud. Using the above HF theory, we show that a 
gas with the same temperature and chemical potential 
as the experimental data is not bose condensed. Thus, a 
mean-field approach doesn't account for our results and 



the quasi-condensate regime is not reached via the usual 
saturation of the excited states. We emphasize that this 
failure of mean-field theory happens in a situation where 
the gas is far from the strong interaction regime, which in 
one-dimensioii (ID), corresponds to the Tonks-Girardeau 
gas limit 0, Q, 0| , and where mean field theory also 
fails. To our knowledge, this is the first demonstration of 
the breakdown of a Hartree-Fock approach in the weakly 
interacting limit. 

We attribute this failure of the mean-field theory to 
the nearly ID character of the gas. It is well known that 
a ID homogeneous ideal Bose gas does not experience 
Bose Einstein condensation in the thermodynamic limit. 
On the other hand, in the presence of repulsive interac- 
tions in the weakly interacting regime, as the linear den- 
sity increases one expects a smooth cross-over from an 
ideal gas regime where q(^)(0 ) ~ 2 to a quasi-condensate 
regime where ff^^HO) l[llEl|l3- The HF approach 
fails to describe this cross-over: as for an ideal gas, the 
thermal fluid does not saturate and g^^^O) = 2 for any 
density. The above results also hold for a ID harmoni- 
cally trapped gas at the thermodynamic limit: no satura- 
tion of excited states occurs for an ideal gas |0| and the 
gas smoothly enters the quasi-condensate regime when 
the peak density increases QJulSj- 

In the experiment presented here, the gas is neither 
purely ID nor at the thermodynamic limit: a few trans- 
verse modes of the trap are populated and a condensation 
phenomenon due to finite size effects might be expected 
|2nf |. However, we will show that, as in the the scenario 
discussed above, the gas undergoes a smooth cross-over 
to the quasi-condensate regime without saturation of the 
excited states. An estimation of the cross-over density 
using a three-dimensional Bogoliubov calculation is in 
agreement with experimental data. 

The experimental setup is the same as [^J. Us- 
ing a Z-shaped wire on an atom chip p^ . we pro- 
duce an anisotropic trap, with a transverse frequency 
of a;^/(27r) = 2.75 kHz and a longitudinal frequency of 
(^2/(277) — 15.7 Hz. By evaporative cooling, we obtain a 
few thousand ^^Rb atoms in the \F — 2,mF = —"2) state 
at a temperature of a few times Tioj^/kB- Current-flow 
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FIG. 1: In situ longitudinal distributions for different temperatures (circles). Long dashed lines: ideal Bose gas profile at 
the shown temperature T and chemical potential /lo obtained from a fit to the wings. Solid lines: profiles obtained in the 
Hartree-Fock approximation for the same T and /lo- Short dashed lines: quasi-condensate profiles with the same peak density 
as the experimental data. Horizontal solid lines: cross-over density estimated using a Bogoliubov calculation (see text). 



deformations inside the micro-wire, located 150 fiui be- 
low the atoms, produce a roughness on the longitudinal 
potential |23. The observed atomic profiles are smooth 
(see Fig.nj, which shows that this roughness is small and 
we neglect it in the following. 

The longitudinal profile of the trapped gases are 
recorded using in situ absorption imaging as in [2lj |. 
The probe beam intensity is about 20% of the satu- 
ration intensity and the number of atoms contained in 
a pixel of a CCD camera is deduced from the formula 
Ns.t — (A^/cr)ln(l2/Ii), where A = 6.0 ^m is the pixel 
size, a the effective cross section, and /i and I2 the probe 
beam intensity respectively with and without atoms. The 
longitudinal profiles are obtained by summing the contri- 
bution of the pixels in the transverse direction. However, 
when the optical density is large and the density varies on 
a scale smaller than the pixel size the above formula un- 
derestimates the real atomic density JlJ. In our case, the 
peak optical density at resonance is about 1.5, and this 
effect cannot be ignored. To circumvent this problem, 
we decrease the absorption cross section by detuning the 
probe laser beam from the F — 2 ^ F' = 3 transition 
by 5 = 9 MHz. We have checked that for larger detun- 
ings, the normalised profile remains identical to within 
5%. For the detuning S, the lens effect due to the real 
part of the atomic refractive index is calculated to be 
small enough so that all the refracted light is collected 
by our optical system and the profile is preserved. 

To get an absolute measurement of the linear density, 
we need the effective absorption cross section a. We find 
a by comparing the total absorption of in situ images 
taken with a detuning d with the total number of atoms 
measured on images taken at resonance after a time of 
flight long enough (5 ms) that the optical density is much 
smaller than 1. In these latter images, the probe beam is 
cr+ polarized and the magnetic field is pointing along the 
probe beam propagation direction so that the absorption 
cross section is 3X^/2tt. We obtain for the in situ images 
taken at a detuning S, an effective absorption cross sec- 
tion (0.24 ± 0.04) X 3A^/27r. For samples as cold as that 
in Fig. n (a), the longitudinal profile is expected to be 
unaffected by the time of flight and we checked that the 
profile is in agreement with that obtained from in situ 



detuned images within 5%. 

Averaging over 30 measured profiles, we obtain a rela- 
tive accuracy of about 5% for the linear density. A sys- 
tematic error of about 20% is possible due to the uncer- 
tainty in the absorption cross section. 

In Fig. n we plot the longitudinal density profiles of 
clouds at thermal equilibrium for different final evapo- 
rating knives obtained from in situ images. We have 
compared the density profiles with the expected quasi- 
condensate density profile with the same peak density. 
This is obtained using the equation of state of the longi- 
tudinally homogeneous gas /i — hLj±\/l + Ana l24l|, and 
the local chemical potential tJi{z) — fio~l/2muj'^z^ , where 
n is the linear atomic density, and hq the chemical poten- 
tial at the center of the cloud. For the two colder clouds 
(graphs (a) and (6)) in Fig.^ we observe a good agree- 
ment between the central part of the experimental curves 
(circles) and the quasi-condensate profile (short dashed 
lines) which indicates that the gas has entered the quasi- 
condensate regime. In (2^, we observed the inhibition of 
density fluctuations expected in this regime. 

For a given atomic density profile, we extract the tem- 
perature and chemical potential of the data from a fit 
to an ideal Bose gas distribution. For cold clouds, the 
ideal Bose gas model should fail in the central part of the 
cloud where interactions are not negligible. We therefore 
fit only the wings of the profile by excluding a number 
of pixels TVox on either side of the center of the profile. 
We fit the longitudinal atomic distribution with only the 
chemical potential /xg as a free parameter for different 
trial temperatures T. As seen in Fig. |21 iJ,Q{T,Nex) is 
approximately linear in Ncx for A^ox > 10. The cloud's 
temperature is the one for which fiQ {T, Ncx) is indepen- 
dent of Ncx- The fluctuations around the straight line 
for Ncx > 10, seen in Fig. |21 is mainly due to the po- 
tential roughness. These fluctuations contribute about 
for 5% to the uncertainty of T and fxo- The determi- 
nation of /io however, is primarily limited by the uncer- 
tainty in the total atom number. For profile (b), we find 
Tfc = (2.75 ± Om)huj^/kB and Mo = (1-65 ± 0.5)huj±. 
For the profile (a), the chemical potential and the tem- 
perature are not very accurate because the wings are too 
small. This profile is not analyzed in the following. 
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FIG. 2: Chemical potential fig, obtained by fitting the wings 
of profile (b) to an ideal Bose gas distribution, as a function 
of the number of excluded pixels A'ex on either side of the 
distribution, for four different trial temperatures. The arrow 
indicates //q.c. ~ hLj±\/l + 4noa, which is the chemical po- 
tential found from the measured peak linear density no. 

Another method to deduce the chemical potential is 
to use the peak atomic density and the formula fi = 
huj±y^l + 4na, assuming the gas is well inside the quasi- 
condensate regime at the center of the cloud. For graph 
(&), we obtain ^'^ ^ = (1.8 ± 0.2) huj±, which is consis- 
tent with the value obtained from fits of the wings of the 
distribution. 

We now compare the experimental density profiles with 
the longitudinal atomic density n{z) given by the Hartree 
Fock theory for the same temperature and chemical po- 
tential. To do the calculation, we also assume the relative 
population of the ground state is negligible (no Bose- 
Einstein condensation) and the quantization of the lon- 
gitudinal eigenstates is irrelevant. In this so called local 
density approximation (LDA), the linear atomic density 
n{z) is identical to the linear atomic density nh(lJ-{z),T) 
of a thermal Bose gas trapped in the radial direction 
and longitudinally untrapped at a chemical potential 
/i(z) — fjLQ — mLo'lz'^ 12. 

To obtain nh{^i,T), we need to compute the self- 
consistent three-dimensional atomic density p which is 
the thermodynamic distribution of independent bosons 
that experience the self consistent HF Hamiltonian|^ Q 

Hhf - + ffkin + Hharm + 2gp{r), (1) 

where iJkin is the transverse kinetic energy term, -ffharm 
the transverse harmonic potential, r is the radial co- 
ordinate and g — ATrh'^a/m where a is the Rb^'' scat- 
tering length. The simplest approach to obtain p{r) 
is to use an iterative method, starting from p(r) = 0. 
For each iteration, we numerically diagonalize the trans- 
verse part of Hhf to deduce the new thermal atomic 
density distribution. When the interactions become too 
strong (linear density larger than 320 atoms per pixels 
for the temperature T^) this algorithm does not con- 
verge. In this case, we use a more time consuming 
method based on a minimization algorithm. We use 
the trial function ptTia.i{r) = ^C2pH2p{r)e^^ /^''owhere 
H2p(r) are Hermite polynomials, and < p < 3. 



We find ro and the four C2p coefficients by minimizing 
i = /o°°'^(piriai('") - Ptriai(r))2rfr//p°°rp2^i^i(r)dr, where 
Atrial ('") ^^"^ thermodynamic equilibrium atomic den- 
sity for iJffF [ptriai(f)]- We find ^ less than 10~^ mean- 
ing that our 5 parameter model describes the trans- 
verse Hartree Fock profile well. The linear density 
rih = J 27rrptriai(''')rf'' is identical to n'^ = J 2Trrp'^^^^^{r)dr 
within 0.5%. In the domain where both methods are 
valid, we also check that they give the same result. 

Figure compares the longitudinal profiles obtained 
with the HF calculation with the experimental data and 
the ideal gas profile for the graphs (5), (c) and (d). The 
HF profile for the hottest cloud (graph (d)) is in agree- 
ment with data and identical within one percent to the 
ideal Bose gas prediction. For the slightly colder cloud of 
graph (c), the HF avoids the divergence in the ideal gas 
model, although it underestimates the peak density by 
approximately 20%. For the even colder cloud of graph 
(6), the discrepancy between the HF profile and the ex- 
perimental data is even larger (35% at the center). 

To validate our Hartree Fock calculations we check a 
posteriori the local density approximation (LDA). The 
LDA is valid if the population difference between adja- 
cent energy states is negligible. This criterion is met if 
the absolute value of pce — Po — ^aip-) is much larger than 
AE, where eo{p) is the ground state energy and AE the 
energy gap between the ground state and the first excited 
longitudinal state. For the temperature of the graph (b), 
as long as fiQ < 2.0huj±, \iJ,es/AE\ is larger than 15 (see 
Fig- El; a-nd for the chemical potential /ig = 1.65Tiu!± de- 
duced from the data \fj,cs/AE\ ~ 25. For such a large 
value of \pcs/AE\, the LDA is expected to be valid. 

For this ratio pcs/AE, we can quantify the error made 
in the density profile due to the LDA. From the HF calcu- 
lation, we obtain the energies i?„(/i,T) of the transverse 
eigenstates. Assuming the transverse motion adiabat- 
ically follows the longitudinal one, we obtain an effec- 
tive longitudinal Hamiltonian with a potential Vn{z) — 
Enipo — rnculz^ /2,T) for each transverse mode. Diago- 
nalization of each effective longitudinal Hamiltonian and 
summation of the resulting thermal profiles gives the ex- 
pected longitudinal density profile. We find agreement 
with the profile obtained from the LDA within 5%. This 
procedure also confirms the LDA for the HF calculations 
corresponding to the graphs (c) and (d). 

In the case of profile (&), where the gas is in the quasi- 
condensate regime at the center, the Hartree-Fock cal- 
culation predicts a population of the ground state A'o ~ 
kBT/{eo - Mo) = 0.0035Ntot- Therefore, the Hartree- 
Fock approach does not predict a saturation of the ex- 
cited states and fails to explain the presence of the quasi- 
condensate at the center of the cloud. The local density 
approximation criterion \pcs\ = \po — £oip)\ >> also 
implies a small relative ground state population. 

The failure of the Hartree-Fock approach for our ex- 
perimental parameters is due to the large density fluc- 
tuations this theory predicts in a dense, nearly ID gas. 
When density fluctuations become too large, pair inter- 
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FIG. 3: Ratio between the effective chemical potential and 
AE, the energy splitting between the ground state and the 
first excited state as a function of the chemical potential for 
the temperature of profile (b). The dashed lines correspond 
to the measured chemical potential. 



actions induce correlations in position between particles 
which are not taken into account in the Hartree-Fock the- 
ory. These correlations reduce the interacting energy by 
decreasing density fluctuations: the gas enters the quasi- 
condensate regime. 

We now estimate the cross-over density rtc.o. at which 
the gas enters the quasi-condensate regime. For this 
purpose, we assume the gas is in the quasi-condensate 
regime and use the Bogoliubov theory to compute den- 
sity fluctuations. We find a posteriori the validity do- 
main of the quasi-condensate regime, which requires 
that density fluctuations 6p be small compared to the 
mean density p. More precisely, we define rico. as 
the density for which the Bogoliubov calculation yields 



JJ{Sp{r))'^/{p{r)nc.o.) (Pr = 1. We indicate this cross- 
over density in Fig^ We find that n^.o. is close to the 
density above which the experimental profiles agree with 
the quasi-condensate profile. 

In conclusion, we have been able to reach a situa- 
tion where a quasi-condensate is experimentally observed 
but a HF approach fails to explain its presence. As 
for purely one dimensional systems, the passage towards 
quasi-condensate in our experiment is a smooth cross- 
over driven by interactions. The profiles that we observe 
require a more involved theory able to interpolate be- 
tween the classical and the quasi-condensate regime. A 
Quantum Monte Carlo calculation th at g ives the exact 
solution of the many body problem [3, Uil should repro- 
duce the experimental data. In fact, since temperatures 
are larger than interaction energy, quantum fluctuations 
of long wavelength excitations should be negligible and 
a simpler classical field calculation should be sufficient 
lis, 26]. Finally, the refined mean field theory proposed 
in [23] , where the coupling constant g is modified to take 
into account correlations between atoms, may explain our 
profiles. 
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